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Particle-like objects are observed to propagate and interact in many spatially extended dynamical 
systems. For one of the simplest classes of such systems, one- dimensional cellular automata, we 
establish a rigorous upper bound on the number of distinct products that these interactions can 
generate. The upper bound is controlled by the structural complexity of the interacting particles — a 
quantity which is defined here and which measures the amount of spatio-temporal information that 
a particle stores. Along the way we establish a number of properties of domains and particles that 
follow from the computational mechanics analysis of cellular automata; thereby elucidating why 
that approach is of general utility. The upper bound is tested against several relatively complex 
domain-particle cellular automata and found to be tight. 
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Persistent, localized, propagating structures — 
particles — have long been observed and constructed in 
cellular automata (CAs) ]l|-p^. A review of the litera- 
ture suggests that particles are widely felt to be some 
of the more interesting phenomena displayed by those 
systems |2^. They are analogous to the "defects" or 
"coherent structures" of pattern formation processes in 
condensed matter physics [^PJ^,|9|, |l6| , |l8| . In fact, in cellu- 
lar automata used to model pattern formation processes, 
the particles model defects and vice versa p^-|27| . A dif- 
ferent analogy to condensed matter physics (specifically, 
hydrodynamics psf) gives them the name "solitons" 
1^,^. They are also known as gliders, glider-like objects, 
or spaceships particularly, but not exclusively, in the 
context of the Game of Life two-dimensional GA [^,^ . 
The name "particle" , while inspired by an analogy to 
field theory in physics, is used here merely for the sake 
of uniform terminology and neutrality of associations.^ 

GA particles, like their physical counterparts, interact 
via collision, and these interactions are well known to 



'Permanent address: Physics Department, University of Wisconsin, Madison, WI 53706. 

^To avoid confusion, we should say that a "particle" in our 
sense is not the same as a par ticle in the sense of interact- 
ing particle systems (IPSs) |j0[|3l| or lattice gases The 
particles of an IPS or the coherent structures that emerge in 
lattice gases may be particles in our sense; we hope to explore 
these and related issues elsewhere. 
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play a crucial role in the dynamics of their underlying cel- 
lular automaton. The construction of computational de- 
vices in CAs, for instance, is almost always accomplished 
through engineering the proper interactions among par- 
ticles ]l| Jl^ , [l3| , ^ -^ . Indeed, it was even at one time 
conjectured by Wolfram |^ that the presence of parti- 
cles in a CA was tantamount to its being computation- 
universal. It is therefore of considerable interest to know 
what interactions a CA's particles may have. Acquiring 
that knowledge is significantly simplified if we can place 
a bound on the number of different interactions between 
any pair of particles. The first successful attempt to do 
so was an expression given in for particles interacting 
on a completely uniform, quiescent background. It has 
been appreciated for some time, however, that many CAs 
display patterned or textured backgrounds — sometimes 
called "domains". Here we substantially generalize the 
original formula to accommodate a large class of domains 
and prove the generalization using elementary automata 
and number theory. 

The outline of the paper is as follows. First, we fix 
basic notation relating to cellular automata, regular lan- 
guages, and finite-state transducers. We then define do- 
mains, in particular periodic domains, and prove some 
basic results about them. We next define particles as a 
particular kind of interface between domains, and define 
interactions between particles. After establishing some 
auxiliary number-theoretic results, we prove the upper 
bound formula that generalizes the main theorem of |^] 
to arbitrary periodic domains. Result in hand, we show 
how it applies to the analysis of several CAs encountered 
in applications, and how it simplifies the analysis of their 
particle dynamics. We close with a summary of the re- 
sults and a list of open questions. An appendix gives the 
details of the proof of an auxiliary result on domains. 

The present work is motivated by and bears on sev- 
eral larger issues. Of particular relevance is the notion of 
an "object" or "coherent structure" that spontaneously 
emerges in the space-time behavior of a process 
The particles analyzed here are arguably one of the sim- 
plest kinds of such emergent structures. Despite this in- 
terest, we do not define "particle" from first principles. 
Like our predecessors we take the existence of particles as 
a given and assume we know how to recognize them in the 
space-time behavior. Nonetheless, the results and their 
proofs do elucidate some of the component concepts that 
we feel will be useful in a theory of emergent structures 
in spatial processes.^ 



II. CELLULAR AUTOMATA, FORMAL 
LANGUAGES, AND TRANSDUCERS 

A cellular automaton (CA) is a discrete dynamical sys- 
tem consisting of a regular lattice of identical cells. At 
each time step t, each of these cells is in one of a number 
k of states E. The state of cell i at time t is denoted 
s\ G E = {0, 1, . . . , fc — 1}. The global state St of a one- 
dimensional CA at time t is the configuration of the entire 



lattice; i.e., sj — (s^ , s 



„i 
t J ■ 



) e , where N is the 



lattice size. One often sees CA phenomenology studied 
where E = {0, 1} and with periodic boundary conditions: 
— ^j^g jnain results reported below do not de- 
pend on these restrictions, however. 

At each next time step t + \, the cells in the lat- 
tice update their states simultaneously according to a 
local update rule cf). This update rule 4> takes as in- 
put the current local neighborhood configuration 77J — 
{sl^^ , . . . ,sl, . . . , sl^^) of cell i and returns the next state 
sl^i; r is the CA's radius. Thus, the CA equations of mo- 
tion are given by 



(1) 



The global update rule $ ; E^ E^ applies in parallel 
(simultaneously) to all cells in the CA lattice; i.e.. 



St+i = $(St) 

= (</.(,7?),</.(^i),...0(7yf-i)) . 



(2) 
(3) 



Binary {k — 2) local state, r — 1 CAs are referred to 
as elementary CAs (EC As) Q . 

An ensemble operator $ can be defined |po| , ^ that 
operates on sets of lattice configurations ilt — {st}: 



such that 



(4) 



(5) 



It is often informative to describe CA configurations 
as one or another type of formal language. A formal lan- 
guage C over the alphabet E is a subset of E* — the set of 
all possible words, or strings, made up of symbols from 
E. A regular language is a formal language whose words 
can be generated or recognized by a device with finite 
memory; sometimes called a finite automaton. Regular 
languages are the simplest class of formal languages in a 
hierarchy (the Chomsky hierarchy) of language classes of 
increasing complexity |Q . 

A deterministic finite automaton (DFA) M is defined 
as a 5-tuple: 



M = {Q,E,<5,go,n 



(6) 



For two approaches to the automatic discovery of particles, 
see I45I and [ftfJl. 
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where Q is a finite set of states, S is an alphabet, qo € Q 
is the initial state, _F C Q is a set of final states, and 
(5:QxE— i-Qisa transition function: 6{q,a) — q' , 
where q,q' € Q and a S S. 

A DFA can be used to read, or scan, words w = 
wi-.-Wl over the alphabet E. Starting in the initial 
state go, the DFA reads the first symbol wi of the word w. 
It then makes a transition to another state q' = S{qo, wi). 
The DFA then reads the next symbol W2 and makes a 
transition to q" = d{q',W2), and so on until all symbols 
in w have been read or until an undefined transition is 
encountered. If, after reading w, the DFA ends in a final 
state q G F, M accepts w; otherwise M rejects it. 

A regular language £ is a formal language for which 
there exists a DFA that accepts all words in C and rejects 
all words not in C. If there is one such DFA, there are 
generally many of them, but there is a unique minimal 
DFA for £, which we write M{C). Similarly, for every 
DFA M there is a corresponding regular language C{M) 
consisting of all and only the words that are accepted by 
M . The regular process languages are a subset of the reg- 
ular languages: those containing all subwords of words in 
the language. All states of the corresponding DFA — the 
process graph — are both initial and accepting states [ pO| . 

Finally, a finite-state transducer (FST) is a finite au- 
tomaton with two kinds of symbol associated with each 
transition: inputs and outputs. Formally, an FST R is 
defined as a 7-tuple: 

R= {Q,^in,^out,5,\,qo,F} , (7) 

where Q, S, qo, and F are as in a DFA, Ei„ is the input 
alphabet, Eout is the output alphabet, and A : Q x — > 
T^out is the observation function: X{q,a) = & where q £ Q, 
a e Sin, and b e Sout- An FST effectively implements a 
mapping fn from one language over Si„ to another lan- 
guage over Eout- In other words, it reads a word w G S*„ 
and transforms it to another word w' e "^out mapping 
each symbol wt £ Sm to a symbol € Sout such that 
w'i — \{q, Wi), where q G Q is the current state of R when 
reading Wi. 

In formal language theory, languages and automata 
play the role of sets and transducers the role of functions. 

III. DOMAINS 

We are now ready to review the computational me- 
chanics analysis of emergent structures in C As |2^j5^ . 

A regular domain A of a CA $ is a process language, 
representing a set of spatial lattice configurations, with 
the following properties: 

1. Temporal invariance (or periodicity): A is mapped 
onto itself by the CA dynamic; i.e., $^A = A for 
some finite p. (Recall that $ takes sets of lattice 
configurations into sets of configurations and that 
a formal language, such as A, is a set of configura- 
tions.) 



2. Spatial homogeneity: The process graph of each 
temporal iterate of A is strongly connected. That 
is, there is a path between every pair of states in 
A'f($'A) for all I. (Recall that M(£) is the minimal 
DFA which recognizes the language C.) 

The set of all domains of a CA <& is denoted A = 
{AO,A\...,A™-i}, where m= [Aj. 

According to the first property — temporal invari- 
ance or periodicity — a particular domain A* consists 
of p temporal phases for some p > 1; i.e.. A* — 
{Ai, A\,..., A;_ J, such that *'Aj = A|^.^,) ^. Here 

p is the temporal periodicity of the domain A'; which we 
denote by T{A^). 

Each of the temporal phases A* of a domain A* is rep- 
resented by a process graph M{A'j) which, according to 
the second property (spatial homogeneity), is strongly 
connected. Each of these process graphs consists of a fi- 
nite number of states. We denote the k^^ state of the j'^ 
phase of A* by A* intentionally suppressing the M{-) 
notation for conciseness. We write the number of states 
in a given phase as S{A'j). 

The process graphs of all temporal phases A^ of all 
domains A* can be connected together and transformed 
into a finite-state transducer, called the domain trans- 
ducer, that reads in a spatial configuration and outputs 
various kinds of information about the sites. (The con- 
struction is given in, for example, |^l|.) Variations on 
this transducer can do useful recognition tasks. For ex- 
ample, all transitions that were in domain Aj's process 
graph are assigned output symbol D, indicating that the 
input symbol being read is "participating" in a domain. 
All other transitions in the transducer indicate deviations 
from the sites being in a domain. They can be assigned 
a unique output ("wall") symbol w G {Wj} that labels 
the kind of domain violation that has occurred. The re- 
sulting domain transducer can now be used to filter CA 
lattice configuration, mapping all domain regularities to 
D and mapping all domain violations to output symbols 
w that indicate domain walls of various kinds. 

We say that a phase of a domain is (spatially) periodic 
when the process graph consists of a periodic chain of 
states, with a single transition between successive states 
in the chain. That is, as one moves from state to state, 
an exactly periodic sequence of states is encountered and 
an exactly periodic sequence of symbols from S is en- 
countered on the transitions. The spatial periodicity of 
a periodic phase is simply ^(A*). We say that a domain 
is periodic when all its phases are periodic. Almost all 
of our attention in the sequel will be confined to peri- 
odic domains, for the following reason. It turns out that 
for such domains all of the spatial periodicities 5 (A*) at 
each temporal phase are equal. Thus, we can speak of the 
spatial periodicity S'(A*) of a periodic domain A*. This 
property, in turn, is central to our proof of the upper 
bound on the number of particle interaction products. 
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Lemma 1 If a domain A* has a periodic phase, then the 
domain is periodic, and the spatial periodicities S{A^^) of 
all its phases A^, j = 0, . . . ,p — 1, are equal. 

Proof. See the Appendix. □ 
Thus, the number of states in the process graph rep- 
resenting a particular temporal phase A^ is the same for 
all j e {1, . . . , r(A*)}, and it is, in fact, 5(A*). 

Finally, there is a larger class of cyclic domains whose 
process graphs consist of a periodic chain of states: as one 
moves from state to state an exactly periodic sequence 
of states is seen. Note that this class includes more than 
periodic domains, which are obviously cyclic. It includes 
domains in which between two successive states in the 
chain there are multiple transitions over E. (See ||5l|| for 
a CA exhibiting two such cyclic domains.) Based on our 
experience we conjecture that Prop. |l| also holds for cyclic 
domains. If this is so, most of the following results, and 
in particular the upper bound theorem, would hold for 
this large class. 

Conjecture 1 For any cyclic domain A\ the spatial pe- 
riodicities S{A'j) of all its phases A^, j = 0, . . . ,p — 1, are 
equal. 



IV. PARTICLES 

When domain violations form a spatially localized (fi- 
nite width), temporally periodic boundary between two 
adjacent domains, they are called particles. 

Definition 1 A particle a is a set {a", a^, . . . , a^^^} of 
finite-width words over S*, called wedges, such that 



*(Aa'A') = Aa(^+i)"^°^PA' 
for some finite p and A and A' G A. 



(8) 



Since a particle is a bounded structure, it does not have 
a spatial periodicity. In the following, the "periodicity of 
a particle" therefore refers to temporal periodicity. 

Since these particles are temporally periodic, we view 
the appearance of wedge as the particle being in it's 
jth phase. The kth symbol in the wedge's word is de- 
noted a^,. The state in which the domain transducer 
finds itself after reading the kth symbol a^, in the wedge 
is denoted q{a'l). 

We now introduce an important but subtle distinc- 
tion. The particle period p referred to above — the surface 
periodicity — is associated with the repetition over time of 
the wedge words as observed in the raw space-time be- 
havior So, Si, S2, . . .. It turns out, as will become clear, 
that particles have an internal periodicity that may be 
some multiple of the surface periodicity p. The internal 
periodicity — the one of actual interest here — though, is 
the periodicity seen by the various phases of the border- 
ing domains. 



Definition 2 A particle a 's intrinsic periodicity P{a) 
is the periodicity of the set of transducer-state sequences 
generated when reading a particle's wedges. For wedge 
= Qfg . . . a^^ the state sequence q^a^) . . . q{a^n) is gen- 
erated in the transducer. We denote this state sequence 
by q{a^). P{a), then, is the number of iterations over 
which the sequence q{a^) reappears. 

Remark 1. P{a) is an integer multiple of a's apparent 
periodicity. 

Remark 2. A simple illustration of the need for intrin- 
sic, as opposed to merely surface, periodicity is provided 
by the 7 particles of EG A 54. See Fig. ||(b) and the 
accompanying text in Sec. VILA, 

After one period P{oi), a particle a will have moved 
a number da of sites in the CA lattice. This shift d^ 
in space after one period is called the particle's displace- 
ment, da is negative for displacements to the left and 
positive for displacements to the right. From the parti- 
cle's periodicity P{a) and displacement da, its average 
velocity is simply Va = da/P{a). 

Note that the above remarks hold whether we are 
looking at the wedges or at the transducer-state labeled 
wedges: one obtains the same velocity. 

The set of all particles a, f3, . . . of a CA $ is denoted 
by P. 

Remark 3. Here we defined temporally periodic par- 
ticles. There are particles in CAs, such as in ECA 18, 
which are temporally aperiodic. In this case, one replaces 
the periodicity condition Eq. (|^) by one using the ensem- 
ble operator; viz., 

*f(AaA') = AaA' . (9) 



A. Structural Complexity of a Particle 

The preceding definitions and discussion suggest that 
one can think of particles as having an internal clock 
or, in the more general case that includes aperiodic par- 
ticles, an internal state, much as the solitary-wave so- 
lutions of continuum envelope equations have internal 
states |2^ . One can ask about how much information a 
particle stores in its states. This is the amount of infor- 
mation that a particle transports across space and time 
and brings to interactions. These considerations lead one 
to a natural measure of the amount of structural com- 
plexity associated with individual particles. 

Definition 3 The structural complexity C{a) of a par- 
ticle a is defined to be 

p-i 

= - E Pr(<7(a^)) log2 Pr(g(a^")) , (10) 
3=0 
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where p is a's period and Pr((7(a^)) is the probabihty of 
a being in phase with the state-sequence q{ce'). 

Remark 1. For the straightforward case of periodic 
particles, in which the wedges and so their associated 
state sequences are equally probable, we have 



C{a) = log2 P{a) 



(11) 



Remark 2. The information available to be processed 
in particle interactions is upper-bounded by the sum of 
the individual particle complexities, since this sum as- 
sumes independence of the particles. As we will see 
shortly, the relative information — that information in one 
particle, conditioned on the other's phase (via the con- 
straints imposed by the mediating domain) and suitably 
averaged — determines the information available for pro- 
cessing by interactions. 



B. Domain Transducer View of Particle Phases 

A particle is bounded on either side by two patches of 
domain. (They could be patches of the same or different 
domains.) Consider what happens to the domain trans- 
ducer as it scans across the part of the lattice containing 
the bounding domains (A* and A* ) and the particle (a). 
It begins by cycling through the states of the process 
graph of a phase (j) of the first bounding domain (A*). 
It then encounters a symbol that does not belong to the 
language of that domain phase, and this then causes a 
transition out of that process graph. Each successive 
symbol of the particle wedge leads to additional transi- 
tions in the transducer. Finally, the transducer reaches 
cells at the beginning of the other bounding domain (A* ) , 
whereupon it begins to follow the process graph of A*, at 
some appropriate phase j' . In this way, a particle wedge 
a' corresponds to a sequence q(a^) of transducer states. 

More formally, the transducer maps a particle wedge 
a\ bordered by A* and A*,, to an ordered n— tuple 
[n = \a^ \ + 2) of states 



Q(a^) = (g(A;.,),g(a^),g(A},,,) 



(12) 



where (/(A* j.) is the transducer state reach on reading 
symbol A^ ^,. Since the transducer-state sequence is de- 
termined by the bounding domain phases and the ac- 
tual wedge , it follows that the mapping from particle 
wedges to state sequences is 1-1. If two particle wedges 
correspond to the same sequence of states, then they are 
the same phase of the same particle, and vice versa. 

This representation of particle phases will prove very 
handy below. 



V. INTERACTIONS 

In many CAs, when two or more particles collide they 
create another set of particles or mutually annihilate. 



Such particle interactions are denoted a -\- [3 — > 7, for 
example. This means that the collision of an a parti- 
cle on the left and a /? particle on the right leads to the 
creation of a 7 particle. Particle annihilation is denoted 
a + (3 ^ %. For completeness we note that there are also 
unstable walls that can spontaneously decay into parti- 
cles. This is denoted a — > /3 -I- 7, for example. 

Often, the actual product of a particle interaction de- 
pends on the phases and (3^ in which the interacting 
particles are at the time of collision. In such a case, there 
can be more than one interaction product for a particular 
collision: e.g., both a -|- /3 — > 7 and a -|- /3 — > can be 
observed. 

The set of a CA's possible particle interactions is de- 
noted I. The complete information about a CA's do- 
mains A, particles P, and particle interactions I can be 
summarized in a particle catalog. The catalog forms a 
high-level description of the CA's dynamics. It is high- 
level in the sense of capturing the dynamics of emergent 
structures. The latter are objects on a more abstract 
level than the original equations of motion and raw (un- 
interpreted) spatial configurations of site values. 



VI. BOUNDING THE NUMBER OF 
INTERACTION PRODUCTS 



Restricting ourselves to particle interactions with just 
two colliding particles — a and /3, say — we now give an 
upper bound on the number ria^fj of possible interaction 
products from a collision between them. (See Fig. ^ for 
the interaction geometry.) In terms of the quantities just 
defined, the upper bound, stated as Thm. ^ below, is: 



< 



P{a)P{(3)Av 
T{Ai)S{A^) 



(13) 



where Av = Va — V/s > and A* is the domain in be- 
tween the two particles before they collide. Note that if 
Av — 0, then ria.p = trivially. 

For simplicity, in the rest of the development we as- 
sume that Aw = Va — Vj} > 0. This simply means that 
particle a lies to the left of (3 and they move closer to 
each other over time, as in Fig. ^ 
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FIG. 1. Interactions between an a and a /3 particle with 
domain A lying between. 



This section proves that Eq. (g_3|) is indeed a proper 
upper bound. The next section gives a number of exam- 
ples, of both simple and complicated CAs, that show the 
bound is and is not attained. These highlight an impor- 
tant distinction between the number of possible interac- 
tions (i.e., what can enter the interaction region) and the 
number of unique interaction products (i.e., what actu- 
ally leaves the interaction region). 

To establish the bound, we collect several intermedi- 
ate facts. The first three lemmas come from elementary 
number theory. Recall that the least common multiple 
lcm(a, b) of two integers a and b is the smallest number c 
that is a multiple of both a and b. Similarly, the greatest 
common divisor gcd(a, b) of two integers a and b is the 
largest number c that divides both a and b. 



Lemma 2 gcd(ca, cb) — c gcd(a, 6), c > 0. 

Proof. See Thm. 2.7 in ||. 

Lemma 3 gcd(a, b) lcm(a, b) = ab. 

Proof. See Thm. 2.8 in ||. 

Lemma 4 lcm(ca, cb) = c lcm(a, b), c > 0. 

Proof. Using Lemmas ^ and ||, it follows that 

lcm(ca, cb) 



□ 



□ 



cacb 



(14) 



icd{ca, cb) 
ab 

= c 

gcd(a, b) 

= c lcm(a, b) . 

□ 

Now we are ready to begin building an analysis of par- 
ticles and particle interactions. 

Lemma 5 The intrinsic periodicity P{a) of a particle a 
is a multiple of the temporal periodicity T{K^) of either 
domain A' for which a is a boundary. That is, 



P{a) = m„,r(A^) , (15) 
for some positive integer rUai that depends on a and A*. 



Proof. At any given time, a configuration containing the 
particle a consists of a patch of the domain A\ a wedge 
belonging to a, and then a patch of A' , in that order 
from left to right. (Or right to left, if that is the cho- 
sen scan direction.) Fix the phase of a to be whatever 
we like — a', say. This determines the phases of A*, 
for the following reason. Recall that, being a phase of 
a particle, a' corresponds to a unique sequence Q{a^) 
of transitions in the domain transducer. That sequence 
starts in a particular domain-phase state A^- ^ and ends in 

another domain-phase state A*, j,, . So, the particle phase 

a' occurs only at those times when A' is in its j*** phase. 
Thus, the temporal periodicity of a must be an integer 
multiple of the temporal periodicity of A* . By symmetry, 
the same is also true for the domain A' to the right of 
the wedge. □ 

Corollary 1 Given that the domain A' is in phase Aj at 
some time step, a particle a forming a boundary of A* 
can only be in a fraction 1/T(A') of its P{a) phases at 
that time. 

Proof. This follows directly from Lemma |^. 

Remark. Here is the first part of the promised re- 
striction on the information in multiple particles. Con- 
sider two particles a and (3, separated by a domain A". 
Naively, we expect a to contain log2 P{oi) bits of infor- 
mation and (3, logj P{0) bits. Given the phase of a, how- 
ever, the phase of hP is fixed, and therefore the number 
of possible phases for fi is reduced by a factor of 1/T(A°). 
Thus the number of bits of information in the a- [3 pair 
is at most 



log2 P{a) -f log2 P(/3) - log2 r(AO) = log, 



P{a)P{P) 
T(A") 



(16) 



The argument works equally well starting from (3. 



Lemma 6 For any two particles a and (3, the quantity 
lcm(P(a), P(/3))Au is a non-negative integer. 

Proof. We know that the quantity is non-negative, since 
the least common multiple always is and Au is so by 
construction. It remains to show that their product 
is an integer. Let ka = lcm(P(a), P(/3))/P(a) and 
kjj — lcm(P(a), P(/3))/P(/3); these are integers. Then 



At; 



2/3 



P{a) 



- kpdp 



lcm(P(a),P(/3)) 



When multiplied by lcm(P(a), P{f3)) this is just kada — 
kpdjs, which is an integer. □ 



6 



Lemma 7 When the distance d between two approach- 
ing particles a and (3, in phases a' and (3^ , respec- 
tively, is increased hy\crsi{P{a),P{f3))/^v sites, the orig- 
inal configuration — distance d and phases and (3^ — 
recurs after lcm{P(a), P{/3)) time steps. 

Proof. From the definition of lcm(a, b) it follows directly 
that lcm(P(a),P(/3)) is a multiple of P(a). Thus, 



(j+lcm(P(a),P(/3))) mod P(q) 



(17) 



and the a particle has returned to its original phase. Ex- 
actly parallel reasoning holds for the /3 particle. So, af- 
ter lcm(P(a), P(/3)) time steps both a and /? are in the 
same phases and again. Furthermore, in the same 
amount of time the distance between the two particles 
has decreased by lcm(pQ,p^)A'y, which is the amount by 
which the original distance d was increased. (By Lemma 
^, that distance is an integer, and so we can meaningfully 
increase the particles' separation by this amount.) Thus, 
after lcm(P(a), P(/3)) time steps the original configura- 
tion is restored. □ 

Lemma 8 // A* is the domain lying between two parti- 
cles a and f3, then the ratio 



lcm(P(a),P(/3))Aw 



(18) 



is an integer. 



Proof. Suppose, without loss of generality, that the par- 
ticles begin in phases a" and f3°, at some substantial 
distance from each other. We know from the previous 
lemma that after a time \cm.{P{a), P{P)) they will have 
returned to those phases and narrowed the distance be- 
tween each other by lcm(P(a), P{(]))Av cells. What the 
lemma asserts is that this displacement is some integer 
multiple of the spatial periodicity of the intervening do- 
main A\ Call the final distance between the particles 
d. Note that the following does not depend on what d 
happens to be. 

Each phase of each particle corresponds to a particu- 
lar sequence of transducer states — those associated with 
reading the particle's wedge for that phase. Reading this 
wedge from left to right (say), we know that Qia'^) must 
end in some phase-state of the domain A'; call it Ag q. 
Similarly, Q(/3°) must begin with a phase-state of A% but, 
since every part of the intervening domain is in the same 
phase, this must be a state of the same phase Aq; call 
it Aq J,. In particular, consistency requires that k be the 
distance between the particles modulo S{A^). But this is 
true both in the final configuration, when the separation 
between the particles is d, and in the initial configuration, 
when it is c? -I- lcm(P(a), P{f3))Av. Therefore 



\cm{P{a),P{(3))Av = d 
\cm{P{a),P{(3))Av = 



(mod SiA')) 
(mod S{A')) 



Thus, lcm(P(Q;), P(/?))Ati is an integer multiple of the 
spatial period ^(A*) of the intervening domain A*. □ 
Remark. It is possible that \cm{P{a), P{P))Av — 0, 
but this does not affect the subsequent argument. Note 
that if this is the case, then, since the least common mul- 
tiple of the periods is at least 1, we have Av = 0. This, 
in turn, implies that the particles do not, in fact, collide 
and interact, and so the number of interaction products 
is simply zero. The formula gives the proper result in 
this case. 

The next result follows easily from Lemmas ^ and ^ 

Lemma 9 // A' is the domain lying between particles a 
and (3, then 

gcd(P(a), P(/?)) = r(A') gcd(m„,, m^,) . (19) 



Proof. We apply Lemma ^: 

gcd(P(a),P(/3)) = gcd(m„,T(A0,m/3,T(A') 
= T(A') gcd{mai,mpi). 



With the above lemmas the following theorem can be 
proved, establishing an upper bound on the number of 
possible particle interaction products. 

Theorem 1 The number Ua^p of products of an interac- 
tion between two approaching particles a and (3 with a 
domain A' lying between is at most 



na,i3 < 



P{a)P{(3)Av 
T(A^)S{A') ■ 



(20) 



Proof. First, we show that this quantity is an integer. 
We use Lemma ^ to note that 

P{a)P{l3)Av _ gcd(P(a),P(/3))lcm(P(a),P(/3))Aw 



T{A')S{A') 



T{A')S{A') 



and then Lemma ^ to find that 

P{a)P{P)Av _ gcd(P(a),P(/3))r 
T{A')S{A^) ~ T(A*) ' 

and finally Lemma ^ to show that 

Pia)P{P)Av ^ r(A')gcd(m,.„m^,)r 
T(A')5'(A^) T(A^) 

= rgcd{mai,mpi) , 



(21) 



(22) 



(23) 



which is an integer. 

Second, assume that, at some initial time t, the two 
particles are in some arbitrary phases a-' and (3-' , re- 
spectively, and that the distance between them is d cells. 
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This configuration gives rise to a particular particle-phase 
combination at the time of collision. Since the global up- 
date function is deterministic, the combination, in turn, 
gives one and only one interaction result. Now, increase 
the distance between the two particles, at time t, by one 
cell, while keeping their phases fixed. This gives rise to a 
different particle-phase combination at the time of colli- 
sion and, thus, possibly to a different interaction result. 
We can repeat this operation of increasing the distance by 
one cell \cm.{P{a), P{f3))Av times. At that point, how- 
ever, we know from Lemma Zjthat after lcm(P(Q!), P{(i)) 
time steps the particles find themselves again in phases 
a' and (3^ at a separation of d. That is, they are in ex- 
actly the original configuration and their interaction will 
therefore also produce the original product, whatever it 
was. 

Starting the two particles in phases a-* and j3-' , the 
particles go through a fraction l/gcd(P(a), P{/3)) of the 
possible P{a)P{fi) phase combinations, over \cm{pa,pp) 
time steps, before they start repeating their phases again. 
So, the operation of increasing the distance between the 
two particles by one cell at a time needs to be repeated 
for gcd(i-'(a), different initial phase combinations. 

This way all possible phase combinations with all possible 
distances (modulo lcm(P(Q!), P(/3))A?;) are encountered. 
Each of these can give rise to a different interaction re- 
sult. 

From this one sees that there are at most 

gcd(P(a),P(/3))lcm(P(a),P(/3))Aw = P(a)P(/3)Au 

(24) 

unique particle-domain-particle configurations. And so, 
there are at most this many different particle interaction 
products, given that <i> is many-to-one. (Restricted to 
the homogeneous, quiescent (A — 0*) domain which has 
T(A) = 1 and S'(A) = 1, this is the result, though not 
the argument, of [^.) 

However, given the phases and (3^ , the distance 
between the two particles cannot always be increased 
by an arbitrary number of cells. Keeping the particle 
phases and (3^ fixed, the amount Ad by which the 
distance between the two particles can be increased or 
decreased is a multiple of the spatial periodicity ^(A*) of 
the intervening domain. The argument for this is similar 
to that in the proof of Lemma ||. Consequently, of the 
lcm(P(a), P{(3))/S.v increases in distance between the two 
particles, only a fraction 1/S'(A*) are actually possible. 

Furthermore, and similarly, not all arbitrary particle- 
phase combinations are allowed. Choosing a phase for 
the a particle subsequently determines the phase A* of 
the domain A* for which a forms one boundary. From 
Corollary d it then follows that only a fraction l/r(A*) 
of the P(/3) phases are possible for the (3 particle which 
forms the other boundary of A'. 

Adjusting the number of possible particle-domain- 
particle configurations that can give rise to different in- 



teraction products according to the above two observa- 
tions results in a total number 



P{a)P{P)/^v 



(25) 



of different particle-phase combinations and distances be- 
tween two particles a and (3. Putting the pieces together, 
then, this number is an upper bound on the number Ua^p 
of different interaction products. □ 

Remark 1. As we shall see in the examples, on the 
one hand, the upper bound is strict, since it is saturated 
by some interactions. On the other hand, there are also 
interactions that do not saturate it. 

Remark 2. We have seen (Corollary |^, Remark) that 
the information in a pair of particles a and /?, separated 
by a patch of domain A*, is at most 



log; 



(26) 



bits. In fact, the theorem implies a stronger restriction. 
The amount of information the interaction carries about 
its inputs is, at most, logj riQ,,/? bits, since there are only 
na,i3 configurations of the particles that can lead to dis- 
tinct outcomes. If the number of outcomes is less than 
"■a,/3j the interaction effectively performs an irreversible 
logical operation on the information contained in the in- 
put particle phases. 



VII. EXAMPLES 



A. ECA 54 and Intrinsic Periodicity 



Figure ^ shows the raw and domain-transducer filtered 
space-time diagrams of ECA 54, starting from a random 
initial configuration. We first review the results of |2|] 
for ECA 54's particle dynamics. 

Figure || shows a space-time patch of ECA 54's dom- 
inant domain A, along with the domain transducer con- 
structed to recognize and filter it out, as was done to 
produce Fig. ||(b). 

Examining Fig. || shows that there are four particles; 
we label these a, (3, 7+, and 7". The first two have 
zero velocity; they are the larger particles seen in Fig. 
H(b). The 7 particles have velocities 1 and —1, respec- 
tively. They are seen in the figure as the diagonally mov- 
ing "light" particles that mediate between the "heavy" a 
and 13 particles. 
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Space 99 

FIG. 3. (a) Space-Time patch of ECA54's primary domain 
A. (b) The transducer that recognizes A and deviations from 
it. After j22). 



The analysis in identified 7 dominant two- and 
three-particle interactions. We now analyze just one: the 
7+ -t- 7~ ^ (3 interaction to illustrate the importance of 
a particle's intrinsic periodicity. 

Naive analysis would simply look at the space-time di- 
agram, either the raw or filtered ones in Fig. ^, and 
conclude that these particles had periodicities ^'(7'*') = 
P{'j~) = 1. Plugging this and the other data — r(A) = 2, 
S'(A) = 4, and Av — 2 — leads to upper bound n^^^ = 
1/4! This is patently wrong; it's not even an integer. 

Figure ^ gives the transducer-filtered space-time dia- 
FIG. 2. (a) Raw space-time diagram and (b) filtered gram for the 7+ and 7~ particles. The domain A is 
space-time diagram of EGA 54 behwior starting from an ar- filtered out, as above. In the filtered diagrams the trans- 
bitrary initial configuration. After [g^. ducer state reached on scanning the particle wedge cells 

is indicated. 
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r-l- 







FIG. 4. The transducer-filtered space-time diagrams for 
tlie 7"*" and 7" particles, (a) The raw space-time patches 
containing the particles, (b) The same patches with the A fil- 
tered out. The cells not in A are denoted in black; those in A 
in white. In the filtered diagrams the transducer state reached 
on scanning the particle wedge cells is indicated. After [ pif . 

From the space-time diagrams of Fig. ||(b) one notes 
that the transducer-state labeled wedges for each particle 
indicate that their intrinsic periodicities are -P(7^) = 2 
and P{-r) = 2. Then, from Thm. we have that 
Tia.p = 1. That is, there is at most one product of these 
particles' interaction. 

Fig. H gives the transducer-filtered space-time diagram 
for the 7+ + 7^ ^ /3 interaction. A complete survey of 
all possible 7+-A-7~ initial particle configurations shows 
that this is the only interaction for these particles. Thus, 
the upper bound is saturated. 



form a class of spatial computations: from all random ini- 
tial configurations, synchronize within a specified number 



of iterations. This CA is 1 



ci of |3q]: a binary, radius-3 



CA. The 128-bit look-up table for (psynci is given in Table 

B 

Here we are only interested in locally analyzing the 
various pairwise particle interactions observed in 4'synci ■ 
It turned out that this CA used a relatively simple set of 
domains, particles, and interactions. Its particle catalog 
is given in Table ||. 

As one example, the two particles a and (3 and the 
intervening domain A have the properties given in Table 
||. From this data, Thm. ^ tells us that there is at most 
one interaction product: 



4-2 



na.j3 



2 • 1 



(27) 



(j) Look-up Table (hexadecimal) 




F8A19CE6B65848EA 




D26CB24AEB51C4A0 


^parent 


CEB2EF28C68D2A04 




E341FAE2E7187AE8 



TABLE I. Lookup tables (in hexadecimal) for (j>synci and 
4>parent- To rccover the 128-bit string giving the CA look-up 
table output bits St+i, expand each hexadecimal digit (the 
first row followed by the second row) to binary. The out- 
put bits st+i are then given in lexicographic order starting 
from the all-Os neighborhood at the leftmost bit in the 128-bit 
string. 




FIG. 5. The transducer-filtered space-time diagrams for 
the 7^ -I- 7~ ^ /3 interaction. After p2|. 



B. An Evolved CA 



The second example for which we test the upper bound 
is a CA that was evolved by a genetic algorithm to per- 
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The single observed interaction between the a and /3 
particles is shown in Fig. ^. As this space-time dia- 
gram shows, the interaction creates another p particle, 
i.e., a + f3 ^ p. An exhaustive survey of the 8 (= 4 x 2) 
possible particle-phase configurations shows that this is 
the only interaction for these two particles. Thus, in this 
case, we see that Thm. |] again gives a tight bound; it 
cannot be reduced. 



Particle Catalog 



Domains A 



Name 


Regular language 


r(A) 


S{A) 


A 


0^0', m* 


2 


1 


Particles P 


Name 


Wall P 


d 


V 


a 


AjA, 4 


-1 


-1/4 




AjAi-j 2 


-1 


-1/2 


7 


A^A,- 8 


-1 


-1/8 


5 


A^A,- 2 








Interactions I 


Type 


Interaction 


Interaction 


React 


a + P~>l3 


J + 13 




React 


S + p^f3 


J + a 


a 


React 


5 + a ^ a 


5 + 7 


ce 


TABLE II 


The particle catalog of <, 


^synci ■ Aj , j 


G {0,1}, 


indicates the two temporal phases of domain A. 






a p 






^ 1 









Time 



P 

99 




Site 99 

FIG. 6. The interaction between an a and a /3 particle in 



C. Another Evolved CA 



The third, more complicated example is also a CA that 
was evolved by a genetic algorithm to synchronize. This 
CA is (pparent of |53|| . It too is a binary radius-3 CA. The 
128-bit look-up table for (jjparent was given in Table ^ 

Here the two particles a and /? and the intervening do- 
main A have the properties given in Table 



[II, Note that 



this is the same domain as in the preceding example. 





^parent 


Particle Properties 




Domain 


T 


S 




A 


2 


1 




Particle 


P 


d 


V 


a 


8 


2 


1/4 


P 


2 


-3 


-3/2 



TABLE III. Properties of two of (^parent's particles. 

From this data, Thm. |l| now says that there are at 
most: 



na.ji 



• 2 • ^ 

4 _i4 



2 • 1 



(28) 



interactions. 

Of these 14 input configurations, it turns out several 
give rise to the same products. From a complete survey 
of a-A-/3 configurations, the result is that there are actu- 
ally only 4 different products from the a + (3 interaction; 
these are: 



a + (3 ^ J 
a + 13 ^2(3 
a + (3 ^ 13 + a 



They are shown in Fig. ^ 

This example serves to highlight the distinction be- 
tween the maximum number of interaction configura- 
tions, as bounded by Thm. ^ and the actual number 
of unique products of the interaction. We shall return to 
this distinction later on. 
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FIG. 7. The four different (out of 14 possible) interaction products for the a + /3 interaction. 



D. EGA 110 



In the next example, we test Thm. |l| on one of the long- 
appreciated "complex" CA, elementary CA 110. As long 
ago as 1986, Wolfram ||l^. Appendix 15] conjectured that 
this rule is able to support universal, Turing-equivalent 
computation (replacing an earlier dictum p^ , p. 31] that 
all elementary CA are "too simple to support universal 
computation" ). While this conjecture initially excited lit- 
tle interest, in the last few years it has won increasing ac- 
ceptance in the CA research community. Though to date 
there is no published proof of universality, there are stud- 
ies of its unusually rich variety of domains and particles, 
one of the most noteworthy of which is Mcintosh's work 
on their tiling and tessellation properties jsj]. Because 
of this CA's behavioral richness, we do not present its 
complete particle catalog and computational-mechanical 
analysis here; rather see ] |55[ ]. Instead, we confine our- 
selves to a single type of reaction where the utility of our 
upper bound theorem is particularly notable. 

We consider one domain, labeled A", and two particles 
that move through it, called /3 and k ]^ . (This /3 particle 
is not to be confused with the (3 of our previous exam- 
ples.) A'' is EC A llO's "true vacuum" : the domain that is 
stable and overwhelmingly the most prominent in space- 
time diagrams generated from random samples of initial 
configurations. It has a temporal period T{hP) = 1, but 
a spatial period ^(A^) = 14. The (3 particle has a pe- 
riod P{P) = 15, during the course of which it moves four 
steps to the left: dp = 4. The k particle, finally, has a 
period P{k) = 42, and moves = 14 steps to the left 
during its cycle. This data gives the f3 particle a velocity 
of vp = —4/15 and the k particle = —1/3. 




Time 



FIG. 8. The particle /3 of EGA 110: The space-time patch 
shows two complete cycles of particle phase. 

Naively, one would expect to have to examine 630 
(= P{(5)P{k) = 15 X 42) different particle-phase con- 
figurations to exhaust all possible interactions. Theorem 
|l|, however, tells us that all but 



(15)(42)(^ - ^) 



(29) 



(14)(1) 

of those initial configurations are redundant. In fact, an 
exhaustive search shows that there are exactly three dis- 
tinct interactions: 

13+ a + SWright , 

P+ K-^ (3 + AWright , 
(3 + K 1-j . 

Here, a, Wright, and rj are additional particles generated 
by EC A 110. These interactions are depicted, respec- 
tively, in Figures [Tl|, 10, and |l2|. 

We should note that the Wright particle is somewhat 
unusual in that several can propagate side by side, or 
even constitute a domain of their own. There are a num- 
ber of such "extensible" particle families in ECA 110 [p5|. 
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Time 




Time 



FIG. 9. The particle k of EGA 110: The space-time dia- 
jram shows one complete cycle of particle phase. 



Time 
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3^right 

Site 



9S 
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FIG. 10. The reaction f3 + k ^ f3 + 4wright in EGA 110. 



FIG. 11. The reaction (3 + k ^ a + Swright in EGA 110. 

We again see, in this complex case, that the bound of 
Thm. |l| is attained. 



VIII. CONCLUSION 
A. Summary 

The original interaction product formula of is lim- 
ited to particles propagating in a completely uniform 
background; i.e., to a domain whose spatial and tem- 
poral periods are both 1. When compared to the rich 
diversity of domains generated by CAs, this is a consid- 
erable restriction, and so the formula does not help in 
analyzing many CAs. We have generalized the original 
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result and along the way established a number of proper- 
ties of domains and particles — structures defined by CA 
computational mechanics. The examples showed that the 
upper bound is tight and that, in complex CAs, particle 
interactions are substantially less complicated than they 
look at first blush. Moreover, in developing the bound 
for complex domains, the analysis elucidated the some- 
what subtle notion of a particle's intrinsic periodicity — a 
property not apparent from the CA's raw space-time be- 
havior: it requires rather an explicit representation of the 
bordering domains' structure. 



Understanding the detailed structure of particles and 
their interactions moves us closer to an engineering disci- 
pline that would tell one how to design CA to perform a 
wide range of spatial computations using various particle 
types, interactions, and geometries. In a complementary 
way, it also brings us closer to scientific methods for an- 
alyzing the intrinsic computation of spatially extended 
systems p5| . 



B. Open Problems 



Time 
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Site 
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FIG. 12. The reaction (3 + k, ^ rj in EGA 110. 



The foregoing analysis merely scratches the surface of 
a detailed analytical approach to CA particle "physics" : 
Each CA update rule specifies a microphysics of local 
(cell-to-cell) space and time interactions for its universe; 
the goal is to discover and analyze those emergent struc- 
tures that control the macroscopic behavior. For now, we 
can only list a few of the open, but seemingly accessible, 
questions our results suggest. 

It would be preferable to directly calculate the number 
of products coming out of the interaction region, rather 
than (as here) the number of distinct particle-domain- 
particle configurations coming into the interaction region. 
We believe this is eminently achievable, given the detailed 
representations of domain and particles that are entailed 
by a computational mechanics analysis of CAs. 

Two very desirable extensions of these results suggest 
themselves. The first is to go from strictly periodic do- 
mains to cyclic (periodic and "chaotic") domains and 
then to general domains. The principle difficulty here is 
that Prop. |l| plays a crucial role in our current proof, but 
we do not yet see how to generalize its proof to chaotic 
(positive entropy density) domains. The second exten- 
sion would be to incorporate aperiodic particles, such as 
the simple one exhibited by EC A 18 [Q. We suspect this 
will prove considerably more difficult than the extension 
to cyclic domains: it is not obvious how to apply notions 
like "particle period" and "velocity" to these defects. A 
third extension, perhaps more tractable than the last, is 
to interactions of more than two particles. The geome- 
try and combinatorics will be more complicated than in 
the two-particle case, but we conjecture that it will be 
possible to establish an upper bound on the number of 
interaction products for n— particle interactions via in- 
duction. 

Does there exist an analogous lower bound on the num- 
ber of interactions? If so, when do the upper and lower 
bounds coincide? 

In solitonic interactions the particle number is pre- 
served |,|,§,|||57|. What are the conditions on the inter- 
action structure that characterize solitonic interactions? 
The class of soliton-like particles studied in ^ possess 
a rich "thermodynamics" closely analogous to ordinary 
thermodynamics, explored in detailed in |^^. Do these 
results generalize to the broader class of domains and 
particles, as the original upper bound of |9| does? 
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While the particle catalog for EGA 110 is not yet prov- 
ably complete, for every known pair of particles the num- 
ber of distinct interaction products is exactly equal to 
the upper bound given by our theorem. This is not gen- 
erally true of most of the CAs we have analyzed and is 
especially suggestive in light of the widely-accepted con- 
jecture that the rule is computation universal. We sus- 
pect that EGA llO's fullness or behavioral flexibility is 
connected to its computational power. (Gf. Remark 2 to 
Thm. |l|) However, we have yet to examine other, com- 
putation universal GA to see whether they, too, saturate 
the bound of our theorem. One approach to this question 
would be to characterize the computational power of sys- 
tems employing different kinds of interactions, as is done 
in I59I ] for computers built from interacting (continuum) 
solitary waves. 
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APPENDIX: PROOF OF LEMMA | 

Lemma 1 If a domain A' has a periodic phase, then the 
domain is periodic, and the spatial periodicities >5'(A* ) of 
all its phases A^, j = 0, . . . ,p — 1, are equal. 

Proof. The proof consists of two parts. First, and most 
importantly, it is proved that the spatial periodicities of 
the temporal phases of a periodic domain A' cannot in- 
crease and that the periodicity of one phase implies the 
periodicity of all its successors. Then it follows straight- 
forwardly that the spatial periodicities have to be equal 
for all temporal phases and that they all must be peri- 
odic. 

Our proof employs the update transducer T^, which is 
simply the FST which scans across a lattice configura- 
tion and outputs the effect of applying the GA update 
rule (j) to it. For reasons of space, we refrain from giving 
full details on this operator — see rather |2^. Here we 
need the following results. If is a binary, radius-r GA, 
the update transducer has 2^^ states, representing the 
2^*" distinct contexts (words of previously read symbols) 
in which scans new sites, and we customarily label 
the states by these context words. The effect of applying 
the GA to a set of lattice configuration represented by 
the DFA M is a new machine, given by T^M — the "di- 
rect product" of the machines M and T^. Once again, 
for reasons of space, we will not explain how this direct 
product works in the general case. We are interested 
merely in the special case where M = Aj, the j*'', peri- 
odic phase of a domain, with spatial period n. The next 



phase of the domain, Aj_^^^, is the composed automaton 
TfjjM , once the latter has been minimized. Before the lat- 
ter step T^M consists of n "copies" of the FST T^, one 
for each of A^ 's n states. There are no transitions within 
a copy. Transitions from copy k to copy k' occur only if 
k' ~ k + 1 (mod n). In total, there are n2^'' states in the 
direct composition. 

T^M is finite and deterministic, but far from mini- 
mal. We are interested in its minimal equivalent machine, 
since that is what we have defined as the representative 
of the next phase of the domain. The key to our proof 
is an unproblematic part of the minimization, namely, 
removing states that have no predecessors (i.e., no in- 
coming transitions) and so are never reached. (Recall 
that, by hypothesis, we are examining successive phases 
of a domain, all represented by strongly connected pro- 
cess graphs.) It can be shown, using the techniques in 
pl[ , that if the transition from state k in A* to state A; -I- 1 
occurs on a (respectively, on a 1), then in the composed 
machine, the transitions from copy k of only go to 
those states in copy k + 1 whose context string ends in a 
(respectively, in a 1). Since states in copy fc -I- 1 can be 
reached only from states in copy k, it follows that half of 
the states in each copy cannot be reached at all, and so 
they can be eliminated without loss. 

Now, this procedure of eliminating states without di- 
rect predecessors in turn leaves some states in copy k-\-2 
without predecessors. So we can re-apply the proce- 
dure, and once again, it will remove half of the remaining 
states. This is because applying it twice is the same as 
removing those states in copy k + 2 for which the last 
two symbols in the context word differ from the symbols 
connecting state k to state k~\-l and state A; + 1 to state 
k + 2 'm the original domain machine Aj . 

What this procedure does is exploit the fact that, in 
a domain, every state is encountered only in a unique 
update-scanning context; we are eliminating combina- 
tions of domain-state and update-transducer-state that 
simply cannot be reached. Observe that we can apply 
this procedure exactly 2r times, since that suffices to es- 
tablish the complete scanning context, and each time we 
do so, we eliminate half the remaining states. We are left 
then with n2^''/2^'' = n states after this process of succes- 
sive halvings. Further observe that, since each state k of 
the original domain machine A* occurs in some scanning 
context, we will never eliminate all the states in copy k. 
Since each of the n copies has at least one state left in it, 
and there are only n states remaining after the halvings 
are done, it follows that each copy contains exactly one 
state, which has one incoming transition, from the previ- 
ous copy, and one outgoing transition, to the next copy. 
The result of eliminating unreachable states, therefore, 
is a machine of n states which is not just deterministic 
but (as we have defined the term) periodic. Note, how- 
ever, that this is not necessarily the minimal machine, 
since we have not gone through a complete minimization 
procedure, merely the easy part of one. A*^j^ thus might 
have fewer than n states, but certainly no more. 
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To sum up: Wc have established that, if A* is a pe- 
riodic domain phase, then A*_|_j^ is also periodic and 
5(A;.+i) < S{A]). Thus, for any t, 5(**Aj.) < S{A^). 
But **Aj) = A\j+t)^„dp and if t = p, we have 

^h+t) mod p = Hi+p) mod p = ^- Putting these together 
we have 

5(A}+i) < S{A)) 5(A}+,) = S{Ai) , (Al) 

for j = 0, 1, ... ,p — 1. This implies that the spatial pe- 
riod is the same, namely n, for all phases of the domain. 
And this proves the proposition when the CA alphabet 
is binary. 

The reader may easily check that a completely paral- 
lel argument holds if the CA alphabet is not binary but 
m— ary, substituting m for 2 and (to — l)/m for 1/2 in 
the appropriate places. We omit it here for reasons of 
space and notational complexity. □ 
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